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The electromagnetic response of graphene, expressed by the dielectric function, and the spectrum 
of collective excitations are studied as a function of wave vector and frequency. Our calculation is 
based on the full band structure, calculated within the tight-binding approximation. As a result, 
we find plasmons whose dispersion is similar to that obtained in the single- valley approximation by 
Dirac fermions. In contrast to the latter, however, we find a stronger damping of the plasmon modes 
Q\ ■ due to inter-band absorption. Our calculation also reveals effects due to deviations from the linear 

Dirac spectrum as we increase the Fermi energy, indicating an anisotropic behavior with respect to 
' the wave vector of the external electromagnetic field. 

<N ' 

5_! ■ PACS numbers: 73.20Mf, 73.22.Lp 



< 

00 



5-1 

CZ5 



O 

o 



I. INTRODUCTION 



Graphene, a single layer of carbon atoms arranged as a honeycomb lattice, is a semimetal with remarkable transport 
properties^. This is due to the band structure of the material which consists of two bands which touch each other at 
' two nodes. The electronic spectrum around these two nodes is linear and can be approximated by Dirac conesi. Many 
unusal transport properties are controlled by the fact that the Fermi energy is at the nodes, where the density of 
states vanishes. However, graphene has been gated such that the Fermi energy can be freely tuned. This has opened 
a wide field for experiments^. In the following we will discuss collective excitation of the 2D electron gas in graphene, 
' caused by an external electromagnetic field. Although this is a standard problem in semiconductor physics, it was 
£h I studied in the case of graphene only in the Dirac approximation around the nodes^i 

I , We consider an electron gas which is subject to an external potential Vi(q, u>). The response of the electron gas 
is a screening potential V s (q, u) which is caused by the rearrangement of the electrons due to the external potential. 
Therefore, the total potential, acting on the electrons, is 



V(c L ,u>)=V i (q,uj)+V s (q,uj) . (1) 



V s can be evaluated self-consistcntly^ and is expressed via the dielectric function e(q, tu). Then the total potential 
£> ' reads^ 

oo : 

V(q,u>) = - rVifaw). (2) 

cn ■ e(q,w) 



The dielectric function is determined by the specific Hamiltonian of the electron gas. In the present case this is the 

. tight-binding Hamiltonian on a honeycomb lattice. 

' The aim of our study is to compare the dielectric function of the full honeycomb lattice with previous studies of 

, the single Dirac cone. In contrast to the latter, we can not rely on a closed expression for the integrals but have to 

^ ' integrate over the Brillouin zone numerically. 
•i— I , 

5-1 1 II. MODEL 

c3 ; 

Using the tight-binding approximation, (quasi-) electrons in graphene are described by the Hamiltonian^ 

H = h\Oi + h 2 <7 2 (3) 

with Fourier components 

3 3 

hi = -i^cos(bj • k), h 2 = -tY^sin(bj • k) . (4) 

i=i 3=1 

CT1.2 arc the Pauli matrices and bi^.3 are the next neighbor vectors on the honeycomb lattice: 



bi = <f(\/3/6, 1), b 2 = d(V3/6,-l), b 3 = d(-l/V3, 0) , 



(5) 
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FIG. 1: Contourplot of the electronic dispersion for different energies. The Brillouin zone is indicated by the hexagon which 
connects the Dirac nodes. 



t is the hopping parameter (ss 2.8eV) and d is the lattice constant (« 1.42 A). The Hamiltonian H satisfies the 
eigenvalue equation H\\d) = £kz|kZ) with eigenvalues E\a = (— l) L E k and the positive square root E k = \fh\ + h\. 
Thus the index I refers either to the conduction band (I = 2) or to the valence band (I = 1). 



III. DIELECTRIC FUNCTION 



The dielectric function can be calculated from the Lindhard formula. Assuming that the wave length of the 
electromagnetic wave is much larger than the lattice spacing, the longitudinal component reads^ 

2?re 2 

e(q,u;) = l x(q,w), (6) 

q 

where 

Xfao,) = jim V = If k J } |(^| e ^|kf)| 2 (7) 

6^o f-f Em - E k >i> + fiu; + ihb 

with k' = k + q. f(E) = l/(e^^ E ~ ,i ' + 1) is the Fermi-Dirac distribution function at inverse temperature (3 = l/fc^T 
with chemical potential \i. In order to determine x(q, uj) we obtain, after some straightforward calculations for Eq. 
([7]), an expression for the polarizability (h~ = 1) as 

X(q, = im 2^ / l K k K k+ q + s - + — — - — , (8) 

6->ajj^ i J BZ \E k - s£/ k+q + uj - id E k -sE k+<l -u + idJ U B z 

where Kk = (hi — ih^jE-^.. flsz is the area of the two-dimensional Brillouin zone Qbz = J BZ cPk. '^' ne i n dex s 
describes intraband (for s = 1) and intcrband scattering (for s = —1). The integral of Eq. (jSJ) is evaluated numerically 
for different values of the chemical potential (i and different directions of the wave vector q, using the limit of zero 
temperature. 



IV. PLASMONS 



Poles in ui of the inverse longitudinal dielectric function l/e(q, uj) for a given wave vector q correspond to collective 
excitations of electrons which are called plasmons. These poles are located either on the real axis or in the complex 
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plane away from the real axis. The latter can be considered as damped plasmons, generated by scattering with 
individual electrons. An imaginary term can appear in the integral x(lj w ) of Eq. ([5]) if the denominator E^—sE^q+tu 
vanishes inside the Brillouin zone. In other words, if (q, ui) is inside the band which is produced by the spectrum of 
the electrons, i.e. where an electronic wave vector k exists that satisfies 

£ k+q -sE w = lj (s = ±1) , (9) 

(cf. Fig. 1), scattering between plasmons and electrons is possible and will lead to damping of plasmons. On the 
other hand, outside the spectrum of electrons (i.e., when there is no electron wave vector k which solves Eq. ((9])) we 
obtain undamped plasmons. 

As an approximation we can expand the electronic dispersion around the nodes in the E = plane. This gives two 
independent Dirac cones (valleys) with linear dispersion E^ ~ 7& around each node. Then the poles of the inverse 
dielectric function can be evaluated exactly^. In this case the plasmon dispersion follows a square root: uip ~ q 1 / 2 . 

Collective excitations, on the other hand, depend on the spectral properties of the electrons. Therefore, deviations 
from Dirac cones can affect them. We found that these deviations lead to a stronger damping of the electrons, 
since electronic excitations require lower energies on the honeycomb lattice in comparison with the linearized (Dirac) 
spectrum. 

In the following we study the imaginary part of the inverse dielectric function^ 

Im ( , 1 , ) = ~ £ /<9 , (10) 
\e(q,w) / e' 2 + e" 2 

where e' (e") is the real (imaginary) part of the dielectric function itself. This quantity becomes a sharp Dirac delta 
function for a pole of l/e(q, ui) on the real axis. If the pole is away from the real axis it becomes a Lorentzian which 
has a width e" , implying that the width of the Lorentzian is a measure for damping by electron scattering. 



V. DISCUSSION 



The polarizability \ of Eq. (JSJ) is evaluated numerically and the corresponding dielectric function is obtained from 
Eq. © for different values of the frequency ui, the wave vector q and chemical potential (Fermi energy) pi. All energies 
are measured in units of the electronic bandwidth (A = 3i), and wave vectors are measured in units of the inverse 
lattice constant. 

The electronic dispersion is plotted for several Fermi energies in Fig. [1] The contours are Fermi surfaces that 
indicate at which wave vectors electronic intraband scattering occurs for a given Fermi energy. For fi = 0.08 we have 
intraband scattering only very close to the Dirac points and the dispersion is Dirac like. As we increase fj, we see the 
effect of "warping" and for fi = 1/3 intraband scattering between different valleys is very likely because we have a 
connected Fermi surface. Finally, for fi = 0.92 we have only a circular Fermi surface around k = with parabolic 
dispersion. This case corresponds with a conventional 2D electron gas. 

The plots in Fig. [2] show the anisotropy of the plasmon dispersion as a consequence of the band structure of 
graphene. Remarkable is the strong deviation from the Dirac case for q x = q y in Fig. (2k- In Fig. [3] it can be seen 
how the plasmon dispersion changes with the Fermi energy \i. The curvature of wp(q) is negative for low energies 
(a square-root behavior close to the Dirac points), it increases with [i and becomes positive for energies higher and 
larger momenta q. This behavior indicates a cross-over from the square-root law of Dirac fermions to the behavior of 
the conventional 2D electron gasi 2 - 



V(4a + v 2 F q) q (q*v F + 4qHi 2 p a + l6k 2 F a 2 ) (v 2 F q + 2a) 
" 4(4a + vj. q )ak F (U) 

with a = 2ne 2 /m, the effective electron mass m, and the Fermi velocity vp. Expansion for small q gives 

3 

w 2 « aq + -v 2 F q 2 , (12) 

which is compared with the findings for graphene at pi = 0.91 (cf. Fig. [3J:). Thus, electrons in graphene behave like 
in a conventional 2D electron gas. This means that if it is possible to vary the Fermi energy between the linear and 
the parabolic regime, the behavior of the fermions in graphene could be 'switched' from relativistic Dirac fermions to 
ordinary electrons. 

In the theory presented above we have made two approximations. First, we have assumed that the wave functions 
4>o(r) of the electrons localized at each atom satisfy the condition 1 0o ( r ) 1 2 ~ <5( r )j i.e. we have ignored the spatial scale 
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ao of <j>o (r). Second, using the Lindhard formulas ©-CO) we have ignored the local field effects^. Both assumptions 
restrict the absolute value of the wave- vector q by the conditions qa® < 1 and qd <C 47r/3, therefore the plasmons 
spectra arc plotted in Figures [2] and [3] at qd < 1.5, which is still within the applicability conditions of the theory. 

In Figure 31 however, we also present our test results for the plasmon spectrum at larger values of q, up to qd ~ 4. 
In spite of at so large q the theory is not quantitatively applicable, one can expect that it provides a reasonable 
qualitative description of the plasmon spectra in graphene. One sees that, while the low-frequency plasmon mode 
with the frequency top ~ Jq gets a large damping and disappears at w/A ~ 0.5, at even higher energies uj/A > 0.8 
the plasmon mode reappears again. The both plasmon branches look like two parts of the same dispersion curve, 
"teared" to the low- and high-frequency pieces by the area of the substantial inter-band damping at 0.5 < uj/A < 0.8. 
This interesting result needs to be studied further within a more advanced theory which takes into account the local 
field effects 1 ^. 



VI. CONCLUSION 

In this paper we have studied the dielectric function of graphene and the plasmon dispersion in tight-binding 
approximation, and compared the results with similar calculations for a single Dirac cone. There are differences due 
to the full energy spectrum of the electrons. It turns out that the damping of the plasmons is stronger, which indicates 
more scattering between plasmons and single electrons on the honeycomb lattice. Moreover, the plasmon dispersion 
is not isotropic. Although the dispersion fits well the sqare-root behavior tup ~ g 1 / 2 at small q values, there are 
substantial deviations for larger values. Finally, a new branch appears in the plasmon dispersion for higher energies 
ui and larger wave vectors q which needs to be subject of further investigations. Our results reflect a crossover from 
Dirac-likc behavior to conventional electron gas behavior by changing the electron density with the help of a gate. 
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FIG. 2: (Color online) Plasmon dispersion for /i = 0.25 and different directions of the q vector a) q x = 0, h)q y — 0, c) q y = q. 
The square-root behavior of the Dirac case is also shown as a dashed curve. 
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FIG. 3: (Color online) Plasmon dispersion for q x = and different values of the Fermi energy /j,. a) /i = 0.33, b) /i = 0.5, c) 
/i = —0.92. In Fig. c) we compare with the plasmon dispersion of the conventional (parabolic) 2D electron gas (dashed curve). 
The curve in c) is the plasmon dispersion of a conventional 2D electron gas of Eq. ([lip, below the dotted line is the intraband 
single particle excitation area of a 2D electron gas. 




FIG. 4: (Color online) Plasmon dispersion for fi = 0.25 and q x = 0. 



